arXiv:1506.01413v2 [cond-mat.str-el] 7 Sep 2015 


Large orbital moment and spin-orbit enabled Mott transition 
in the Ising Fe honeycomb lattice of BaFe 2 (P 04)2 

Young-Joon Song^, Kwan-Woo Lee^’^Q and Warren E. Pickett^ 

^Department of Applied Physics, Graduate School, Korea University, Sejong 339-700, Korea 
^ Department of Display and Semiconductor Physics, Korea University, Sejong 339-700, Korea 
^Department of Physics, University of California, Davis, CA 95616, USA 
(Dated: September 8, 2015) 

BaFe 2 (P 04)2 is an unusual Ising insulating ferromagnet based on the Fe^^ spin S = 2 ion, the 
susceptibility of which suggests a large orbital component to the Fe local moment. We apply density 
functional theory based methods to obtain a microscopic picture of the competing interactions and 
the critical role of spin-orbit coupling (SOC) in this honeycomb lattice system. The low-temperature 
ferromagnetic phase displays a half-semimetallic Dirac point pinning the Fermi level and preventing 
gap opening before consideration of SOC, presenting a case in which correlation effects modeled 
by a repulsive Hubbard U fail to open a gap. Simultaneous inclusion of both correlation and SOC 
drives a large orbital moment in excess of 0.7 p,B (essentially L = 1) for spin aligned along the 
c axis, with a gap comparable with the inferred experimental value. The large orbital moment 
accounts for the large Ising anisotropy, in spite of the small magnitude of the SOC strength on the 
3d (Fe) ion. Ultimately, the Mott-Hubbard gap is enabled by degeneracy lifting by SOC and the 
large Fe moments, rather than by standard Hubbard interactions alone. We suggest that competing 
orbital occupations are responsible for the structural transitions involved in the observed re-entrant 
rhombohedral-triclinic-rhombohedral sequence. 

PACS numbers: 71.20.Be, 71.30.-l-h, 75.25.Dk, 71.27.-|-a 


I. INTRODUCTION 

In the last few decades, the two-dimensional (2D) 
honeycomb lattice has attracted interest due to several 
exotic physical phenomena related with superconduc¬ 
tivity, magnetism, and topological phases, and to the 



FIG. 1: (Color online) (a) The rhombohedral structure, which 
contains PO 4 tetrahedra and edge-sharing FeOe octahedra, of 
BaFe 2 (P 04)2 in its low and high T regimes, (b) Top view of 
the two-dimensional Fe honeycomb lattice in the a — b plane. 
The solid line indicates the unit cell, containing two Fe ions 
(brown circles). The oxygen ions (small red circles) in the 
upper layer are denoted by a dot. 
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Dirac point 2D material graphene. Recently, discover¬ 
ies of tcmological insulators, originally on the honeycomb 
lattice, [y have stimulated increased research on proper¬ 
ties of systems with honeycomb lattices. M 

Transition-metal systems are good examples to in¬ 
vestigate the interplay among lattice, spin, orbital de¬ 
grees of freedom, and correlation effects, since variation 
of 3c? —>■ 4(i —>■ 5c? leads to weakening strength of corre¬ 
lation, but strengthening of spin-orbit coupling (SOC). 
For 5c?^ systems, a large SOC results in a Jeff = I 
Dirac (or relativistic) Mott insulator in Ba2Na0s06[l, @ 
in which the orbital moment plays a pivotal role. Sim¬ 
ilarly, a Jeff = 5 Mott transition has been proposed 
in 5c?® Sr2lr04, i.e., one hole in the Og manifold. Q 
4c?^ systems usually have a reduced SOC and some en¬ 
hanced correlation strengths. In 3d systems, a much 
stronger correlation strength leads to a Mott transition, 
which has been intensively discussed, 0 while SOC is 
usually minor. Here, we focus on the recently synthe¬ 
sized BaFe2(P04)2 compound with a honeycomb lat¬ 
tice of high spin d® Fe^'*' ions, leading to an effectively 
isolated minority-spin d^ configuration due to a large 
exchange splitting of the high spin Fe ion. Large 3d 
orbital moments have been encountered recently[9| in 
low-symmetry environments but gaining understanding 
of their origin is not simple. 

Mentre and coworkers synthesized a rare example 
of an insulating 2D Ising ferromagnetic (FM) oxide 
BaFe2(P04)2, confirmed by the observed critical expo¬ 
nents in good agreement with the theoretical Ising val¬ 
ues, with the Curie temperature Tc = 65 K and negligi¬ 
ble hysteresis. This system also shows an unusual 

re-entrant structural transition sequence: a re-entrant 



TABLE L Optimized internal parameters and structure infor¬ 
mation with the experimental lattice parameters of a=4.869 
A and c=23.230 A at low r.[^ In the representation of a 
hexagonal lattice, the Ba ions are at 3a (0,0,0). The Fe, Of, 
and P atoms sit at 6c (0,0,z), while 02 lies at 18/ {x,y,z). 
The 02 ions are corner shared by the PO 4 tetrahedra and the 
FeOe octahedra, whereas the 01 ions are at one of the ver¬ 
tices of the PO 4 tetrahedra. The bond length d and the bond 
angles are given in units of angstrom and degree, respectively. 


Fe P 

internal 2 ; 2 

parameter 0.1680 0.5710 

01 02 

z X y z 

0.6371 0.3494 0.0235 0.8840 

d(Fe-0)(x3) 


2.042, 2.050 

d(P-O) 


1.536, 1.598(x3) 

ZO-Fe-0(x3) 


86.6, 88.6, 93.2,91.8 

ZO-P-O (x3) 


107.2, 111.7 


rhombohedral (i?3) —>■ triclinic (Pi at 140 K) —>■ rhombo- 
hedral ferromagnetic (FM) (P3 at 70 K). This rare reen¬ 
trant sequence was suggested to arise from com pet ition 
between Jahn-teller distortion and magnetism. |10l ini 
The observed Curie-Weiss moment, 6.16 is sub¬ 
stantially enhanced from the spin-only value of 4.9 ^b 
for a S = 2 system. Additionally, the saturated mo¬ 
ment is about 5 ^lb from powder neutron-diffraction 
studies. ©[HI These values are close to those of the case 
of S' = 2 , L = 1 , thus requiring a large orbital moment on 
the 3d ion. The latter corresponds to the possible largest 
orbital moment Mb = 1 /tb for the system. Optical 
spectroscopy measurements show semiconducting behav¬ 
ior with an estimated energy gap of 1.5 eV. [3 

Although the generalized gradient approximation 
(GGA) plus on-site Coulomb repulsion U calculations 
by Mentre and coworkers provided some preliminary 
results, 0 available information about the electronic 
structures is still limited. Here we report more extensive 
density functional theory (DFT) based studies, including 
correlation and SOC, to uncover the origin of the strong 
magnetocrystalline anisotropy, Mott transition, and the 
large orbital moment in BaFe 2 (P 04 ) 2 . 


II. CRYSTAL STRUCTURE AND METHODS 

Figure[T]shows the rhombohedral structure of the space 
group R3 (no. 148), which is observed in both low and 
high T regimes. In this structure, P ions sit in the cen¬ 
ter of the Fe honeycomb lattice, but above and below 
the layer along the c direction. The P ions are edge 
shared with three oxygen ions of the FeOe octahedra, and 
form PO 4 tetrahedra. The Fe^"*" layers are isolated by 
the intervening insulating (P 04 )^“ tetrahedra and 
ions. In the low T phase, the interlayer distance is 7.71 
A along the c axis, while the Fe-Fe distance is 2.81 A. 
Both structurally and electronically (below), the com¬ 


pound is strongly 2D in nature. 

Using the GGA exchange-correlation functional, the 
internal parameters (iM [l3 | were optimized at the experi¬ 
ment lattice parameters for the low T phase (see Table [I]), 
since no experimental data are available. The optimized 
positions are very similar to the experimental ones at 
high T, as expected since the difference in volumes be¬ 
tween the low and high T phases is only several tenths 
of a percent. PEI The FeOe octahedra are somewhat 
more regular in the low T phase, while the PO 4 tetrahe¬ 
dra are close to an ideal tetrahedron in the high T phase. 

Our calculations were performed in a rhombohedral 
unit cell with GGA exchange-correlation functional plus 
U and SOG, which are implemented in the accurate 
all-electron full-potential code Wien2k,[i^ which uses 
the tetrahedron method for Brillouin-zone integrations. 
The fully localized limit flavor of the GGA-I-C/ functional 
was used for this strong Fe moment compound. On-site 
Goulomb repulsion U = 3 — 5 eV and the Hund’s ex¬ 
change parameter J = 0.70 eV were employed, as widely 
used for Fe^"*" systems. 00 The basis size was deter¬ 
mined from the smallest atomic sphere by RmtKmax = 7 
and augmented plane-wave radii of Ba 2.30, Fe 2.00, P 
1.61, and O 1.29, in atomic units. The Brillouin zone 
was sampled up to a fc mesh of 21 x 21 x 21 to check con¬ 
vergence carefully in this system, which is nonmetallic at 
almost all levels of the theory. 


III. RESULTS 

A. Underlying electronic structure 

Before discussing the low-temperature FM phase re¬ 
sults, we address the nonmagnetic state, which repre¬ 
sents the underlying electronic structure. The enlarged 
band structure in the energy range of Fe d orbitals and 
the corresponding densities of states (DOSs) are given in 
Fig. [21 In the region displayed there are ten d bands, 
since the unit cell contains two Fe ions. From the tiny 
oxygen character in the DOS around Ep^ we conclude 
that direct Fe Og ~ Og hopping is substantial due to the 
edge-sharing octahedral structure, and this interaction 
leads to a bonding-antibonding splitting ~I eV of the 
fully filled t 2 g manifold. A gap of 0.25 eV opens between 
the antibonding t^g state and the unfilled eg manifold. 
This feature of isolated bands around Ep in a hexagonal 
lattice system is similar to that of the hole-doped super¬ 
conductor Lii_a;Nb 02 . 00 In both bonding and anti¬ 
bonding manifolds, Dirac points appear at the K point 
by symmetry, as observed in several honeycomb lattices. 

An anomalous feature is that the top of the t^g man¬ 
ifold (bottom of the gap) is flat along the M — K line, 
hence along the entire edge of the Brillouin zone (which 
consists of alternating M — K and K — M lines). Gon- 
sidering also the dispersionless character along the c axis 
due to the large interlayer distance, the electronic struc¬ 
ture becomes an unusual quasi-one-dimensional one, as 
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FIG. 2: (Color online) Nonmagnetic GGA band structure 
(left) and orbital projected densities of states (right) of the 
rhombohedral structure in the energy region containing Fe d 
orbitals, not that the 01 contribution is negligible and invis¬ 
ible here. The occupation confirms the d® Fe^^ configuration 
and the density of states just below the gap shows a quasi- 
one-dimensional character due to the flat band that extends 
from M to K, hence all around the edge of the Brillouin zone. 
The bottom of the gap is denoted by the horizontal dashed 
line. The fc-point labels are for the conventional hexagonal 
lattice, viz. graphene. 

evident in the DOS near Ep. This one-dimensionality, 
however, is not uniaxial but instead radial, perpendic¬ 
ular to the edge of the Brillouin zone. A similar one- 
dimensionality was found as well in the hexagonal lattice 
CuAlOs compound. 0 


B. Ferromagnetic electronic structure 

1. Uncorrelated electronic structure and magnetism 

Now we focus on the FM state, which is the ground 
state at the low T regime. The magnetic energy is large, 
with the FM state being 1.2 eV/Fe lower in energy than 
for the nonmagnetic state. This very large magnetization 
energy guarantees that the local moment is robust even 
in the high T disordered regime. Based on the Stoner 
picture (which may be accurate only for small moments) 
where this energy is |/M|, Ms = d/is leads to a broad 
first estimate of the Stoner I of 0.3 eV. A more realistic 
means of estimating I is from the band exchange splitting 
IMs = Aex = 3 eV of the GGA bands, giving I = 0.75 
eV, more in line with accepted values for the magnetic 
3d ions. 

The energies of FM alignment and in-plane antiferro¬ 
magnetic (AFM) alignment have been compared, varying 
the value of the on-site repulsion U. For U equal to 3 eV 
or less, FM alignment is favored. For 17=4-5 eV, however, 
AFM alignment is favored. In this range the orbital occu¬ 
pations for FM alignment do not change, whereas there 



FIG. 3: (Color online) (a) FM band structures and DOSs 
within GGA in the range of Fe d bands, indicating a half- 
semimetallic character in the minority bands (red dashed 
lines). A Dirac point appears almost exactly at Ep at the 
point K. (b) Overlaid FM GGA (red dashed lines) and 
GGA-bSOC (green solid lines) band structures and total DOS 
enlarged near Ep, in a ±0.3-eV region containing only some 
of the minority t 2 g manifold. The Dirac point is split by about 
30 meV by SOG. Note that a tiny hole pocket compensating 
a tiny electron pocket at the K point appears in a region that 
is not visible along the lines shown. 


is some small evolution of orbital occupations for AFM 
order. The orbital moment is the same large value in all 
cases. Evidently this question of magnetic alignment is 
near a delicate balance, so we have not pursued this ques¬ 
tion of trying to obtain the nearest-neighbor exchange 
coupling either in plane (which this difference reflects), 
or out of plane, where coupling should be much smaller 
and also involves the difficult question of frustrated cou¬ 
pling. We proceed with analysis using the observed FM 
alignment. 

The left side of Fig. [3Ka) displays the GGA FM band 
structure, showing a half-semimetallic character with a 
minority-spin Dirac point (or Weyl node)[23 lying right 
below Ep by 30 meV at the point K. The corresponding 
DOS, given in the right side of Fig. UDa), shows a i- 
hlled minority t 2 g manifold, confirming the high spin S 
= 2 Fe^"*" configuration. The t 2 g manifold shows a large 
exchange splitting of 3 eV. All majority states {t 2 g and 
Bg) are filled, hence | spin, and orbitally inert. Hence¬ 
forward we focus on the single minority electron, whose 
spin half reduces the total spin to S' = 2 and introduces 
the question of orbital occupation. 

With the local X, Y, Z coordinates of the FeOe octa¬ 
hedron, the t 2 g orbitals with global z axis along c can be 
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FIG. 4: (Color online) Band structures near Ef: (a)-(c) 
within GGA+U, and (d)-(f) from GGA+17+SOC, for U = 
1 — 4 eV and J = 0.7 eV. The spin-up and -down bands in 
(a)-(c) are denoted by black solid and red dashed lines, re¬ 
spectively, while the bands including SOC (d)-(f) are shown 
by green solid lines. In (d) the scale is enlarged to allow the 
gap opening at the Dirac point to be seen, (g) Charge density 
plot of the minority d orbitals, shown at the isovalue of 0.09 
e/A®. The c axis is along the vertical direction. This shape is 
nearly identical for all calculations performed here, (h) total 
DOS at [/ = 4 eV in GGA+SOC-l-I/, showing a 1.2-eV energy 
gap between the corresponding up and down bands. 


described in a three-fold symmetry adapted manner by 

= -^{CndxY + ^ladyZ + ^^dzx), ( 1 ) 

where the phase factor is = exp{i^^), m = is the 
projection of the orbital moment, and the superscript is 
an exponent. Due to the small trigonal distortion, the 
Lz = 0 orbital (po {dz^ in shape) has a somewhat higher 
site energy than the doublet 4>±i. The corresponding 
DOS near Ep has some similarity to a particle-hole asym¬ 
metric version of graphene, [2l| due to a Dirac point at 
K that lies very near the Fermi level. The strong mag¬ 
netism completely changes the Fe ion configuration from 
to 



U(eV) 


FIG. 5: (Color online) Variation of the Fe orbital Ml moment, 
and energy gap Eg, varying U in GGA-|-[/-|-SOC from 0 to 6 
eV, with the spin directed along the c axis. From U = 4 eV, 
the gap is between the Filed spin-up and the unfilled spin- 
down channels. The critical value of U is in the 2.5-3-eV 
range. 


2. Dirac point: Semimetal-insulator transition 

We now address the observed insulating character of 
BaFe 2 (P 04 ) 2 , using the GGA-bD approach for correlated 
insulators. Applying a small U = 1 eV with J = 0.5 — 
1 eV to the Fe ions within GGA-I-C/, the Dirac point 
degeneracy pins Ep and no gap opens, as shown in Fig. 
SKa). The Dirac cone is isotropic, as in graphene. As 
U is increased, the Dirac point degeneracy continues to 
pin Ef. Applying U up to 7 eV, beyond the limit of a 
proper value for Fe ions, although higher energy bands 
shift this system remains a half-semimetal without gap, 
as illustrated in Fig. |3Kb) and (c). 

As mentioned in the Introduction, experiment shows 
semiconductin g b ehavior and implies a large orbital 
moment Mp , E [HI which are not explained by 
GGA-bG. We carried out calculations including SOC, 
both GGA-bSOC, and GGA+C/-bSOC. The GGA-bSOC 
band structure, overlapping the GGA one, is shown in 
the left panel of Fig. EDb). SOC leads to a band splitting 
of 40 meV near Ep, visible in this plot only at the F and 
K points. The linear Dirac bands are also split by 30 
meV. This small effect is consistent with the small SOC 
strength in Sd ions, and spin mixing is also degraded by 
the 2.5-eV exchange splitting. 

Another interesting feature occurs in the parabolic 
bands touching Ep at the F point. Inclusion of SOC 
splits the bands, degenerate at the F points, and leaves 
only one maximum ai Ep, which is more isotropic than 
in the case excluding SOC. As a result, as shown in the 
DOS of the right side of Fig. [SKb), a van Hove singu¬ 
larity (vHs) appears very near Ep, in addition to two 
additional vHs at -0.15 and 0.25 eV. 

Simultaneous inclusion of correlation and SOC effects 
splits the Dirac point degeneracy, thereby opening a gap 
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for U as small as 1 eV. At small U this gap is a spin-orbit 
gap, not a Mott gap. Although the strength of SOC is 
small in Sd systems, the symmetry lowering due to SOC 
becomes crucial for o pen ing a gap in this system, as oc¬ 
curs also in BaCrOa.j^ The resulting band structures 
in GGA-bt/-bSOC are given in Fig. |4i;d)-(f). As dis¬ 
played in Fig. [5l the energy gap of 1.2 eV for [/ = 4 
eV is close to the experimentally estimated value. The 
gap is indirect, between K and F, although due to the 
small dispersion the direct gap is not much larger. There 
has been some uncertainty in determining the optical gap 
spectroscopically. [I^ 


3. Degeneracy lifting and large orbital moment 

GGA-bSOC gives a negligible magnetocrystalline 
anisotropy, with an orbital moment of Ml ~ O.OS/xs/Fe 
both for in-plane and c-axis spin alignment. Varying U in 
GGA-I-G-bSOC with the spin along c. Ml of Fe rapidly 
increases in the U=2-3-eV range, then saturates near 
O.T/xs at [7 = 4 eV, as displayed in Fig. [S] This value 
excludes contributions from the tails of the 3d orbitals 
extending out of the Fe sphere, so the atomic value will 
be somewhat larger. (This large orbital moment arises 
for antiferromagnetic alignment as well.) As a critical 
value Uc of 2.5-3 eV, the opening of a Mott gap begins. 

Consistent with the large exchange splitting A^x = 3 
eV, the total spin moment Ms = d/is/Fe remains un¬ 
changed, since mixing of the spin-split orbitals by small 
SOC is minor. At U = 4 eV, the spin direction along the 
c direction has much lower energy around 400 meV/f.u. 
than for spin in the d — b plane, explaining the observed 
Ising behavior. The state obtained for the spin direction 
in the d — b plane (not pictured) remains semimetallic 
(Dirac point) even in the large U regime, showing a very 
small Ml ot '^0.06 fiB- This striking anisotropy conhrms 
that SOC drives the semimetal-insulator transition, and 
that depends strongly on spin direction. This type of 
Ising magnet, i.e., with the spin direction perpendicular 
to the layers, has been proposed as a crucial ingredient 
for emergence of a quantum anomalous Hall stateQ in 
ferromagnets. 

We now discuss the origin of the large orbital moment. 
The SOC operator is 

Hso=iS-L = + \{S+L_ + S.L+)]. (2) 

The majority shell is filled, spin split from the active 
minority orbital by A^x, and inert, because spin mixing 
is reduced by ^/Aex Thus the spin and orbital 

operators are those of the minority channel. With the 
majority channel frozen out, the only effect is from the 
diagonal term SzL^ with Sz = — 5 . For the magnetiza¬ 
tion along the c axis, the orbital degeneracy is broken by 
— ^Lz for the two orbital projections. With a bandwidth 
of 0.5 eV, this small splitting induces an orbital moment. 


but more importantly it specifies the orbital to be occu¬ 
pied in the Mott insulating phase. For the magnetization 
in plane, the Tz = ±1 degeneracy remains unbroken, no 
orbital moment arises, and correlation is ineffective in 
opening a gap, as discussed above. 

As can be seen in Fig. HKg), the spin-density plot of the 
filled minority t 2 g manifold shows the <f)^i character char¬ 
acteristic of an Lz = 1 orbital. When symmetry allows 
it, the two orbitals (j)±i are equally filled, the semi-Dirac 
degeneracy remains at AT, and so no gap emerges. It is 
SOC that breaks the degeneracy, which though small, en¬ 
ables the correlation effects to operate. Upon increasing 
U, the single minority electron density begins to transfer 
from the if-i orbital to the (j)+i orbital. 


IV. DISCUSSION AND SUMMARY 

Using the DFT-based calculations, we have investi¬ 
gated the Ising insulating Fe honeycomb lattice system 
BaFe 2 (P 04 ) 2 , which at the GGA level results in a half- 
semimetallic electronic structure with a fermionic Dirac 
point lying nearly precisely at Ep; this degeneracy is 
split by inclusion of SOC. Applying both U and SOC 
suppresses the direct t 2 g — t 2 g interaction, and leads to a 
Mott transition of uncommon origin, becoming enabled 
by the symmetry lowering due to the (small in magni¬ 
tude) SOC. In the insulating phase the magnetic phase 
is S' = 2, L = 1. The atomic-like value of L suggests 
classifying BaFe 2 (P 04)2 as a Je// = 3 Dirac Mott insu¬ 
lator. These DFT-based studies have explained the large 
magnetic anisotropy, the very large orbital moment for 
a 3d ion, and the unusual manner in which spin-orbit 
coupling enables the Mott gap to open in BaFe 2 (P 04 ) 2 . 
Even the magnitude of the gap is given reasonably by 
our methods. 

The re-entrant structural transitions mentioned in the 
Introduction needs comment. The Fe honeycomb lat¬ 
tice studied here bears a great deal of similarity to the 
(lll)-orient ed p erovskite bilayers that are recently be¬ 
ing studied,!^, a difference being that in the lat¬ 
ter case the metal-oxygen octahedra are corner sharing 
rather than edge sharing. In those systems structural 
symmetry breaking is predicted in some cases, using the 
same methods as used here. That symmetry breaking 
involves an unusual type of orbital-occupation competi¬ 
tion. The high-symmetry structure (threefold, as here) 
supports complex-valued, orbital moment containing or¬ 
bitals, but breaking to low symmetry (PI space group in 
that case) is connected with occupation of t 2 g type real 
orbitals adapted to the local pseudocubic symmetry. 

The same orbital-occupation competition can be ex¬ 
pected in BaFe 2 (P 04 ) 2 . At high temperature where the 
Fe spins are disordered (primarily up and down due to the 
Ising nature) the (space- and time-averaged) symmetry is 
the high threefold symmetry, and the symmetry adapted 
(j)+i orbital is occupied and shows up in the Gurie-Weiss 
moment. The “ground state” (with spins still disordered) 
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may be one in which a real valued t 2 g orbital is ordered 
instead, with accompanying symmetry lowering, and the 
transition toward such a lower-symmetry structure oc¬ 
curs at 140 K. However, upon ordering ferromagnetically, 
the strains that accompany occupation of cubic orbitals 
(viz. Jahn-Teller strains) may be incompatible with long- 
range order. Such spin-lattice competition could then fa¬ 
vor a return to the symmetric structure at or near the 
Curie temperature, as observed. Further consideration 
of this mechanism will be left for future work. 

A Dirac point has played an important role in our 
discussion of how the Mott gap in BaFe 2 (P 04)2 arises 
only when spin-orbit coupling is included. Few in¬ 
stances of a Dirac point pinned to the Fermi level in 
a ferromagnet have been reported, which produces a 
half-semimetallic system. One such instance is in the 
Cr 02 /Ti 02 multilayer [ 2 ^, where it is actually a semi- 
Dirac semi-Weyl point. Another instance is that of (111)- 
oriented perovskite superlattices, [U, [ 2 ^ where Dirac 
point degeneracies are only lifted by symmetry-breaking 
interactions. Related band points arise in the hexago¬ 
nal W lattice on the Cl-Si surface, [41 and are also of 
interest in triangular lattice systems. 27| In some such 
systems, Chern bands with a Chern number ±2 are 


proposed. [ 2 ^[ 2 ^ In addition, the three-dimensional Weyl 
semimetals NbP and TaAs, which are distinguished by 
lack of inversion symmetry but preserving time reversal- 
symmetry, have been very much of interest, since these 
compounds show an extremely large and unsaturated 
magnetoresistivity even for very high magnetic field. [ 2 ^ 
[3l| BaFe 2 (P 04)2 bears some similarity to the above- 
mentioned compounds, and has the potential for addi¬ 
tional interesting properties. 
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